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We study flavour susceptibilities in the A/" = 4 SU{N) super Yang-Mills plasma coupled to two 
quark flavours at strong coupling and large A'^ by means of its gravity dual. The off-diagonal 
susceptibility is 1/A'^-suppressed and we compute it as a one-loop effect on the gravity side. Contrary 
to naive extrapolation from perturbative results, it attains a finite value in the infinite-coupling limit. 
Moreover, its parametric form is independent of whether or not meson bound states exist in the 
plasma. We conclude that caution must be exercised when drawing conclusions about the QCD 
plasma from lattice calculations of quark susceptibilities. 



1. Introduction. At a temperature ~ 170 MeV, 
Quantum Chromo dynamics (QCD) undergoes a cross 
over into a deconfined phase [l[ , the so-called quark-gluon 
plasma (QGP). Asymptotic freedom ensures that a de- 
scription in terms of weakly coupled quarks and gluons 
is applicable at T 3> Tc. An outstanding question is 
whether a similar description in terms of weakly cou- 
pled quasiparticles is possible in the vicinity of Tc- Al- 
though this regime can be rigorously studied in lattice- 
regularized QCD, in this formulation the nature of the 
underlying degrees of freedom is not directly apparent 
but instead must be inferred from the behaviour of com- 
putable quantities. It is thus important to check the 
validity of these inferences in gauge theory plasmas in 
which the underlying degrees of freedom are understood 
both at weak and at strong coupling. 

Computable quantities in finite-temperature lattice 
QCD consist mainly of thermodynamic potentials and 
charge susceptibilities. Unfortunately, the former are 
rather insensitive to the nature of the effective degrees 
of freedom. This is clearly illustrated by the four- 
dimensional N = A super Yang-Mills (SYM) plasma, 
which has proven to be a good toy model for the QGP just 
above Tc (see e.g. [H and references therein). At weak 
coupling the SYM plasma is well described as a gas of 
gluons, adjoint scalars and fermions. At strong coupling 
no quasiparticles exist and the effective description in- 
volves black hole in five-dimensional anti-dc Sitter (AdS) 
space 3 . Despite this dramatic difference, in the infinite- 
coupling limit the pressure saturates at 75% of the free 
value [4|. Remarkably, the pressure of QCD in the region 
of interest deviates from the free limit by a comparable 
fraction [l[ . In addition, perturbative calculations gener- 
ically suffer from problems when extrapolated to such 
low temperatures and the theoretical uncertainties be- 
come comparable to the deviations above @ . 

In contrast to thermodynamic potentials, the flavor 
structure of the QGP has been regarded as a better dis- 
criminant for the presence of quasiparticles. Quark sus- 



ceptibilities are defined as 
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with Z the QCD partition function and fii = Hu.d 
u and d chemical potentials. Early lattice calculations 
of diagonal quark susceptibilities [6[ indicated a rapid 
convergence to the free-theory limit above . This were 
contrasted with a relatively larger deviation (by a factor 
of 2) between the free and the strongly coupled limits 
of the A/" = 4 SYM plasma. However, more recent lattice 
results show a slower approach to the free limit, with 
sizable deviations up to temperatures of 2 Tc. As before, 
in this regime the uncertainty arising from perturbative 
calculations [l^ becomes comparable to the difference 
between the free and the strongly coupled limits [1] . 

Similarly, off-diagonal quark susceptibilities have also 
been proposed [ll| as a sensitive probe of the degrees of 
freedom of the plasma, based on the assumption that a 
dominance of quark-like quasiparticles must lead to un- 
corrclated flavors and hence to small off-diagonal suscep- 
tibilities. In perturbative SU{N) QCD, the off-diagonal 
susceptibility is stron gly suppressed with respect to the 
diagonal one because [10| 
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Xuu 



with g the coupling constant and X = g^N the 't Hooft 
coupling. In contrast, lattice results have shown that this 
ratio is sizable up to temperatures as high as 2 Tc. Since 
this cannot be fully accounted for by the perturbative 
result it has been taken as an indication of the pres- 
ence of bound states in the plasma [l2j (see also (l3|). 
In contrast to the previous cases, this inference cannot 
be tested in strongly coupled A/" = 4 SYM coupled to 
dynamical quarks because the off-diagonal quark suscep- 
tibility in this theory has not been calculated. The reason 
is that, as can be seen from the off-diagonal suscepti- 
bility is 1/A^-suppressed in the large- limit, and hence 
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it corresponds to a quantum effect on tfie dual gravity 
side. The purpose of this Letter is to perform this calcu- 
lation and thus shed light on the interpretations of the 
lattice results described above. We will focus on deter- 
mining the qualitative form of the result as well as its 
leading parametric dependence. 

2. Gravity calculation. The M = ^ SYM plasma cou- 
pled to two quark flavours {u and d) of equal masses is 
dual to two coincident D7-brane probes [IJ] in AdS^ x 
in the presence of a black hole horizon [3|. At temper- 
atures T < M the branes sit entirely outside the hori- 
zon in the so-called 'Minkowski phase' (MP) [isj. The 
spectrum in the flavour sector consists of mesons with 
masses M„„3 ~ M and parametrically heavier quarks 
with masses ~ VXM On the gravity side 

the quarks correspond to long strings stretching between 
the branes and the horizon, and the mesons to normal- 
mode fluctuations of the branes. Both sets of excita- 
tions are good quasiparticles since their widths vanish 
in the limit iV, A — >■ oo. In contrast, at temperatures 
T > M the branes pierce the horizon in the so-called 
'black hole phase' (BHP). In this phase no quasiparticles 
exist. Meson-like excitations are now broad resonances 
whose frequencies posses large imaginary parts, and cor- 
respond on the gravity side to quasinormal modes on the 



brane [17[ . Quarks are effectively massless in this phase, 



but the widths of quark-like excitations are comparable 
to their energies. We will determine the susceptibilities 
in both phases. 

The embedding of the D7-branes in AdS^ x in the 
presence of the black hole depends only on the dimen- 
sionless ratio M/T. Following we write the in- 

duced metric on the branes (in Euclidean signature) as 
ds^ = L'^ds^ig), with 

ds'^ig) = grrdr'^ + gxxdxf + grrdr^ + gnndfll . (3) 

The gauge theory coordinates x = {r, x'} and the AdS 
radial coordinate r are dimensionless. Dimensions can 
be restored through the replacements x — > ttTx, r 
t/ttL^T. The last term in ([3]) is the metric on the S'^ 
wrapped by the branes. The metric functions depend 
on M/T and r, which we choose so that r G (0, cxj) in 
both phases. In the BHP g-rr vanishes at the horizon, 
i.e. grr = at r = 0. 

The full action describing the non-Abelian U{2) gauge 
theory on the branes is not known. However, for our 
purposes it suffices to consider the Yang-Mills action [l^ 



S 



drd''x^TY{F''''F^,} 



(4) 



where the trace is over flavour indices and includes 
gQQ. The effective coupling constant is = Att'^/N [i8| . 
so the large- TV limit corresponds to the classical limit. 

We introduce chemical potentials by turning on a back- 
ground value for the gauge field on the branes 



Br 



-i (mbI + f^y) H{r) , 



(5) 



where = {^u + Md)/27rr and /^i = (/x„ - /.td)/27rT 
are the (dimensionless) baryon and isospin chemical po- 
tentials, repectively, and is the two-dimensional Pauli 
matrix. The function H{r) is obtained by solving the 
classical equations, which yields H{r) = 1 in the MP 
and 



J r 

in the BHP, where 
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is a dimensionless function that depends only on the cm- 
bedding of the branes. Substituting ([5]) into the action 
(jl]) we obtain (the potential-dependent part of) the free 
energy, defined as F = TS^ — —TlogZ. The result is 
Fo = in the MP and 

^o = -^/(f)^(M^+M^) (8) 



in the BHP, where the superscript '0' indicates that this is 
the classical result. Using ([T|) we get the susceptibilities: 

, [MP] (9) 

[BHP] (10) 
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The diagonal values reproduce the original calculation 
[20t . We will come back to these results in Sec. 3. Suffice 
it to note here that the off-diagonal susceptibilities van- 
ish at leading order in in both phases, as expected. We 
will compute the leading-order correction in the 1 /N ex- 
pansion. The corrections that are sensitive to the chem- 
ical potentials arise from quantum fluctuations of the 
brane flelds around the classical solution ([5|). For con- 
creteness we will focus on vector fluctuations, i.e. we set 



(11) 



The matrix-valued field a can be decomposed into modes 
with different charges under the U (2) fiavour group as 

(12) 



The Lagrangian C for these excitations is obtained by 
substituting eqn. pT|) into the action ^ . To capture the 
leading 1 /N corection we expand the lagrangian to lead- 
ing order in e, i.e. to quadratic order in a. Due to the 
commutation relations among Pauli matrices, the only 
components that couple to the background field (HI are 
For simplicity, henceforth we restrict ourselves to 
the two transverse excitations that satisfy k^a]''^ = 
in Fourier space. Since the Euclidean time is periodic 
under r — > r -I- tt, each of these can be expanded as 
a° = e~2^*a°, so we finally arrive at 
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(13) 
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where we have introduced ai = {a\ +iaj ) /^/2. For 
Minkowski embeddings, for which H{r) = 1, this is the 
Lagrangian of a scalar with a chemical potential Hi . 

To determine the partition function, we will express a 
general fluctuation a in term of a basis of (radial) func- 
tions obtained by solving the eigen- value problem 



1 



\ (14) 
Normalizability requires the eigen-functions 0£„ to vanish 
as r — >■ oo. As r — > regularity requires — > in the 
MP [11], and r^l^l in the BHP [2lj. In addition, 

these functions obey the orthonormality conditions 
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In this basis, the /ii-independent part of the Lagrangian 
is diagonal, which yields a Gaussian partition function 
for ^1 = given by 

= n / ^^*endcin exp [-c|„C£„ (fc^ + k^„)] . (16) 

ink •' 

The full partition function is obtained by averaging the 
/ii-dependcnt part over the partition function ()16p : 



exp 



= Z(0) X 

^ ' inm 



(17) 
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where V is the 3-dimcnsional volume and 



^inl = / dr^gg^^g-^H{r)4>}^<j>,^ , (18) 
rfrVff .9"^9""i?(r)'(/<L</'^m ■ (19) 



(2) 



Since we are only interested in susceptibilities, we can 
consider infinitesimal values of p^i/T and determine the 
partition function via perturbation theory. Since at van- 
ishing /i, the plasma is isospin neutral, the leading-order 
perturbation in /ii/T is given by 



(51ogZ = /ify 
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This expression only depends on and on the mode 
spectrum, which is determined by the branes' embed- 
ding, i.e. by M and T alone. In particular, it is both A- 
and A'^- independent, irrespectively of the phase. Taking 
derivatives with respect to fi^ and /i„ we conclude that 
the contribution of the transverse vector modes takes the 
form 

xL = Xdd = -Xld ^Hf)T\ [MP& BHP] (21) 



where /i is a coupling- and A'-indepcndent function, and 
the superscript '1' indicates that this is the one- loop re- 
sult. We will come back to it in Sec. 3. 

The interpretation of eqn. (PO)) is obscured by the fact 
that the Euclidean solutions ([Tl)) do not correspond to 
the physical excitations of the plasma, i.e. to mesons in 
the MP or to wide resonances (quasinormal modes) in the 
BHP. Both of these are defined in Lorentzian signature, 
and they are characterized by dispersion relations 0Jn{k) 
and the associated radial wave-functions ipn (k). The con- 
nection can be made clear by noticing that the Euclidean 
differential equation ([T4|) and the corresponding bound- 
ary conditions are obtained from their Lorentzian coun- 
terparts upon the analytic continuation 



±iKen , w 



-2i£. 



(22) 



This implies that, viewed as functions of complex mo- 
mentum, the dispersion relations obey a;„(±i«;£„) = 
— 2z£, and similarly that (^„(±iK£„) = (/)£„. FoUowig (2l| 
we rewrite eqn. (j20p by noting that its poles coincide with 
those of 

d^k v-V- 



= M^^/t^I.!. (23) 

4e^^^'('=) jv£^(fc) ^ 32£2 e^^i^(^)A^£2(fc)jvW (fc) 
|iw„(fc) - 2£|' ^ |zw„(A:) - 2^1' \iLo„-,{k) - 2e\\ 

where Pn\k) and Pmn{k) are polynomial functions and 



dr^gg--g-''H{r)^l{k,r)ipr.{k,r) , (24) 
dr^gg^^g--H{rf^l{k,r)^^{k,r) (25) 



have no poles in the complex plane. 

The coefficients (f24| -(|25 |) now have a simple interpre- 
tation. The radial wave functions of the physical exci- 
tations satisfy the orthonormality conditions (different 
from (HH)) 



dry/g g'''' g'-'''iPnik,r)(f^{k,r) = Sn 



(26) 



Therefore the functions Nnm{k) are the matrix elements 
of the operator H(rY in this basis. It is clear that these 
mixing coefficients determine the (second order) pertur- 
bation of the physical modes in the presence of the ex- 
ternal field. 

In the MP the function H{r) is unity and hence 

Nnm{k) — 5nm. Also, in this case it is easy to see [ll| 
that the polynomials P^^'^{k) vanish. Since the frequency 
of the modes does not depend on the Matsubara frequen- 
cies 2£, eqn. ([23]) coincides with that of a charged bosonic 
field. Performing the sum over £ we then arrive at the 
one-loop contribution from transverse vector mesons 

— ^ ^, [MP] (27) 



TJ (271)3(1 
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where we have restored physical units and one factor of 
2 arises from the two transverse polarizations. We con- 
clude that in the MP there is no mode mixing, i.e. mesons 
contribute to the partition function independently of one 
another. In contrast, in the BHP H{r) is not constant 
and the mixing coefficients Nnm are not diagonal. As a 
consequence, the contributions to the one-loop suscepti- 
bilities of the broad resonances present in the plasma mix 
with one another. 

3. Discussion. The A/" = 4 SYM plasma provides a 
rich model in which several expectations about the be- 
haviour of quark susceptibilities at strong coupling can 
be tested. One important conclusion is that, contrary 
to what naive extrapolation of the perturbative formula 
^ might suggest, the off-diagonal susceptibility attains 
a finite value in the limit A — > oo, as illustrated by 
eqn. (PT|) . On the gravity side this follows simply from 
the overall factor 1/e^ in the action (|4]), together with 
the quadratic nature of the leading-order Lagrangian for 
the fluctuations. Since these features are common to all 
brane fluctuations (bosonic or fermionic) , Xud will remain 
A-independent once the contributions of all the fluctu- 
ations are included. This conclusion applies regardless 
of the nature of the effective degrees of freedom in the 
plasma. Indeed, the MP contains quark-like and meson- 
like quasiparticles, whereas no quasiparticle excitations 
of either type exist in the BHP. Despite this dramatic 
difference, the off-diagonal susceptibilities are parametri- 
cally the same in the two phases. In particular, the BHP 
illustrates the fact that non-vanishing off-diagonal sus- 
ceptibilities do not imply the existence of narrow meson- 
like states. This limited sensitivity suggests that caution 
should be exercised when extracting conclusions about 
the value of the coupling constant or about the nature 
of the underlying degrees of freedom in the QGP from 
lattice calculations of the off-diagonal susceptibilities. 

The fact that the off-diagonal susceptibility is of or- 
der 7V° in both phases is due to the fact it arises from 
contributions of meson-like states in both cases. In the 
MP these are narrow quasiparticles and eqn. pT]) shows 
that each of them contributes independently, as expected 
In the BHP these are broad resonances and their 
contributions mix with one another. 

Diagonal susceptibilities are of order 7V° in the MP 
and of order TV in the BHP. In the MP this is due to the 
fact that quark-like excitations are parametrically heav- 
ier than the mesons by a factor of V^, so their contribu- 
tion to the susceptibilities is not visible at leading order. 
In contrast, in the BHP quarks are effectively massless 
and they produce an order-A contribution to the diago- 
nal susceptibilities, eqn. (fTO|). despite the fact that they 
do not give rise to narrow quasiparticles. In the BHP 
eqn. ([21]) provides the oder-one meson correction to the 
leading quark contribution. 

The change in scaling of the diagonal susceptibilities 



produces a change in the ratio Xud/Xuu from order one 
in the MP to order 1 /N in the BHP. We stress that this 
strong decrease in the ratio is not due to a reduction of 
cross-flavour correlations across the transition but rather 
to the appearance of additional, flavour-charged degrees 
of freedom (the quarks) in the BHP. 
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